Phase transitions in diluted negative- weight percolation models 
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We investigate the geometric properties of loops on two-dimensional lattice graphs, where edge 
weights are drawn from a distribution that allows for positive and negative weights. We are interested 
in the appearance of spanning loops of total negative weight. The resulting percolation problem is 
fundamentally different from conventional percolation, as we have seen in a previous study of this 
model for the undiluted case. 

Here, we investigate how the percolation transition is affected by additional dilution. We consider 
two types of dilution: either a certain fraction of edges exhibit zero weight, or a fraction of edges is 
even absent. We study these systems numerically using exact combinatorial optimization techniques 
based on suitable transformations of the graphs and applying matching algorithms. We perform 
a finite-size scaling analysis to obtain the phase diagram and determine the critical properties 
of the phase boundary. We find that the first type of dilution does not change the universality 
class compared to the undiluted case whereas the second type of dilution leads to a change of the 
universality class. 

PACS numbers: 64.60.ah, 75.40.Mg, 02.60.Pn, 68.35.Rh 



I. INTRODUCTION 

Percolation is one of the most fundamental problems 
in statistical mechanics [1, Q . Its ideas apply to physi- 
cal disciplines as different as solid-state physics and cos- 
mology and it has the ability to describe many phase 
transitions in a completely geometric and probabilistic 
context. The pivotal question in percolation theory is 
that of connectivity. Consider for example random-bond 
percolation where one studies a lattice that is randomly 
filled with a certain fraction of edges. Clusters of con- 
nected sites are then analyzed regarding their geometric 
properties. Depending on the fraction of edges on the 
lattice, the geometric properties of the clusters change, 
leading from a phase with rather small and disconnected 
clusters to a phase where there is basically one large 
cluster dominating the lattice. Therein, the appearance 
of an infinite, i.e. percolating, cluster is described by a 
second order phase transition. In the past decades, a 
large number of percolation problems have been investi- 
gated in various contexts through numerical simulations. 
Among those are problems where the objects of interest 
are string-like, rather than clusters as in the standard 
case explained above. The geometric properties of string- 
like objects have been studied and quantified in diverse 
contexts @, i, 0, @, EE El EE El • Recently [H we have 
investigated negative weight percolation (NWP), a prob- 
lem in similar vein but with subtle differences compared 
to other string- like percolation problems. In NWP, we 
consider a regular lattice with periodic boundary condi- 
tions where all adjacent sites are joined by undirected 
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edges. Weights are assigned to the edges, representing 
quenched random variables drawn from a distribution 
that allows for edge weights of cither sign. The details 
of the weight distribution is therein controlled by a dis- 
order parameter. For a given realization of the disorder, 
we then ask for a configuration of loops, i.e. closed poly- 
gons on the lattice, such that the sum of the edge weights 
that make up the loops is minimal. The problem of find- 
ing these loops can be cast into a minimum-weight-path 
(MWP) problem as will be explained below in more de- 
tail. This MWP is also of relevance for the problem of 
finding domain-wall excitations in 2d Ising spin glasses 
EH . One feature of this particular MWP problem is that 
there is no modularity, meaning that a path obeying the 
optimization constraint can have subpaths that do not 
meet the optimization criterion. We further impose the 
constraint that the loops are not allowed to intersect, as 
a result there is no definition of clusters in NWP. 

Previously we investigated NWP on two-dimensional 
lattices. We found a critical value of the disorder param- 
eter above which percolating loops emerge in the limit 
of large system sizes. We further determined a whole 
set of critical exponents that characterize the underlying 
disorder-induced and geometric phase transition. These 
exponents where clearly different from those describing 
other percolation transitions. Moreover, the critical ex- 
ponents where found to be universal in 2c?, meaning that 
they do not depend on the details of the disorder distri- 
bution or the lattice geometry. 

Albeit there are essential differences between NWP 
and usual bond percolation described above, the features 
of the latter are important if we ask for the effect of dilu- 
tion on the characteristics of NWP. Here, we investigate 
the effect of dilution on the critical properties of NWP. 
Therefore we use observables from percolation theory and 
a finite-size-scaling (FSS) analysis to probe critical lines 



2 



(a) (b) 




percolating loops 



FIG. 1: Samples of L = 64 loop configurations for the two 
different types of disorder investigated. Percolating loops ap- 
pear as dashed lines and are marked by arrows, (a) p — O.f, 
pi = 0.3940, where p\ is the fraction of edges with weight equal 
to zero. A loop can include those edges without increasing its 
configurational energy. Edges with zero weight that are in- 
cluded in loops are colored grey, (b) p — 1.0, pu = 0.4998, 
where pu is the fraction of missing edges on the lattice. Those 
edges that comprise the cluster of connected edges on which 
the percolating loop is located appear as thin grey lines. 

in the disorder-dilution plane for two different types of 
dilution described below. Further, we check whether the 
critical exponents change along the critical lines. We 
therefore analyze the critical exponent v that describes 
the divergence of a typical length scale in NWP as we 
approach the critical point. For one distinguished point 
in the disorder-dilution plane, a more complete analysis 
involving more critical exponents is given as well. 

The paper is organized as follows. In section [Til we 
introduce the model in more detail and we outline the 
algorithm we used to obtain the loop configurations. In 
section lliTl we list the results of our numerical simulations 
followed by a summary and conclusions in section ITVl 

II. MODEL AND ALGORITHM 

In the framework of this paper, we consider 2d square 
lattices G— (V,£) with side length L and fully periodic 
boundary conditions (BCs), i.e. N = L x L sites i e V 
and a maximal number of 2N (in case of no dilution) 
edges {i,j}&£ joining adjacent sites i,j€V. We further 
assign a weight /cost u>ij to each edge contained in £, 
representing quenched random variables that introduce 
disorder to the lattice. One realization of the undiluted 
disorder therein consists of a fraction p of edge weights 
drawn from a Gaussian distribution with zero mean and 
unit width, together with edges with unit weight (fraction 
1 — p) . This weight distribution explicitly allows for loops 
C, i.e. closed polygons on the lattice, with a negative 
total weight toe = j}ec ^ or anv nonzero value of 
the disorder parameter p, a sufficiently large lattice will 
exhibit at least small loops with a weight smaller than 



zero. If the disorder parameter is sufficiently large, even 
large loops with negative weight will appear, which span 
the system. Hence, the system is percolating in this case. 
Here, as an additional source of disorder, we consider two 
different types of dilution, affecting the lattice edges: 

• Type I: The disorder is characterized by a fraction 
pi of edge weights wy = 0. This does not mean that 
the respective bonds are absent in terms of bond 
percolation. Instead, a loop can include these edges 
without an increase of the configurational energy, 
see figure Ufa); 

• Type II: The lattice is diluted by a fraction pu 
of absent bonds. Depending on the value of pu, 
this induces a usual random bond percolation pro- 
cess that will affect the negative weight percolation 
(NWP) of loops, see figure[]Jb). 

Given G together with a realization of the disorder, 
we determine a set C of loops such that the configura- 
tional energy, defined as the sum of all the loop-weights 
E = X)z;gc CJ ' c ' l& minimized. As further optimization 
constraint, the loops are not allowed to intersect and 
generally, the weight of an individual loop is smaller 
than zero. Note that C may also be empty. Clearly, 
the configurational energy E is the quantity subject to 
optimization and the result of the optimization proce- 
dure is a set of loops C, obtained using an appropriate 
transformation of the original graph as detailed in [i~2| . 
So as to identify the edges that constitute the loops for 
a particular instance of the disorder, we need to obtain a 
minimum- weig htcd perfect matching (MWPM) [H, [H 
on the transformed graph. 

Here, we give a brief description of the algorithmic 
procedure that yields a minimum weight set of loops for 
a given realization of the disorder. Figure O illustrates 
the 3 basic steps, detailed below: 

(1) each edge, joining adjacent sites on the original 
graph G, is replaced by a path of 3 edges. Therefore, 
2 "additional" sites have to be introduced for each edge 
in £. Therein, one of the two edges connecting an ad- 
ditional site to an original site gets the same weight as 
the corresponding edge in G. The remaining two edges 
get zero weight. The original sites i G V are then "dupli- 
cated", i.e. i — > zi,?2) along with all their incident edges 
and the corresponding weights. For each of these pairs 
of duplicated sites, one additional edge {ii,^} with zero 
weight is added that connects the two sites i\ and ii. 
The resulting auxiliary graph Ga = (Va, £ a) is shown in 
figure [DJb), where additional sites appear as squares and 
duplicated sites as circles. FigureHJb) also illustrates the 
weight assignment on the transformed graph Ga- Note 
that while the original graph (figure [U^a)) is symmetric, 
the transformed graph (figure [^b)) is not. This is due 
to the details of the mapping procedure and the particu- 
lar weight assignment we have chosen. A more extensive 
description of the mapping can be found in [To| . 

(2) a MWPM on the auxiliary graph is determined 
via exact combinatorial optimization algorithms [l|. A 
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FIG. 2: Illustration of the algorithmic procedure, (a) origi- 
nal lattice G with edge weights, (b) auxiliary graph Ga with 
proper weight assignment: black edges carry the same weight 
as the respective edge in the original graph and grey edges 
carry zero weight, (c) minimum-weight perfect matching M: 
bold edges are matched and dashed edges are unmatched, 
(d) loop configuration (bold edges) that corresponds to the 
MWPM depicted in (c). 

MWPM is a minimum- weighted subset M of £a, such 
that each site contained in Va is met by precisely one 
edge in M. This is illustrated in figure [He), where the 
solid edges represent M for the given weight assignment. 
The dashed edges are not matched. Due to construction, 
the auxiliary graph consists of an even number of sites 
and since there are no isolated sites, it is guaranteed that 
a perfect matching exists. Even if theres no edge on G, 
i.e. pu = 1 in terms of the dilution introduced above. In 
that extreme case, Ga consists solely of the "duplicated" 
sites that are joined pairwise by zero-weighted edges. A 
perfect matching trivially exists and the corresponding 
MWPM has weight zero. 

(3) finally it is possible to find a relation between the 
matched edges M on Ga and a configuration of negative- 
weighted loops C on G by tracing the steps of the trans- 
formation (1) back. In this regard note that each edge 
contained in M that connects an additional site (square) 
to a duplicated site (circle) corresponds to an edge on G 
that is part of a loop, see figure [^d). More precisely, 
there are always two such edges in M that correspond 
to one loop segment on G. All the edges in M that con- 
nect like sites (i.e. duplicated-duplicated, or additional- 
additional) carry zero weight and do not contribute to a 
loop on G. Once the set C of loops is found, a depth-first 
search [12, [HI can be used to explore C and to determine 



the geometric properties of the individual loops. For the 
weight assignment in figure Ufa), there is only one loop 
with weight ujc — — 2 and length I = 8. 

Note that the result of the calculation is a collection C 
of loops such that the total loop weight, and consequently 
the configurational energy E, is minimized. Hence, one 
obtains a global collective optimum of the system. Ob- 
viously, all loops that contribute to C possess a nega- 
tive weight. Regarding the weight assignment in step (1) 
there are different possibilities that all result in equiv- 
alent sets of matched edges on the transformed lattice, 
corresponding to the minimum- weight collection of loops 
on the original lattice. Some of these weight assignments 
lead to a more symmetric transformed graph, see e.g. 
[l2| . However, this is only a technical issue that does 
not affect the resulting loop configuration. Albeit the 
transformed graph is not symmetric, the resulting graph 
(figure [5£d)) is again symmetric. 

In the following we will use the procedure outlined 
above to investigate the effect of dilution on the critical 
properties of NWP. 



III. RESULTS 

So as to clarify the effect of dilution on the critical 
properties of NWP, we first need to locate the criti- 
cal lines in the p — pi/n-planes that separate the per- 
colating from the nonpercolating domains. Therefore, 
we use observables from percolation theory and a finite- 
size-scaling analysis. Besides the weight ujc and length 
£ = Y^{ij}ec 1 °f an individual loop, we determine its 
linear extensions by projecting it onto the perpendicular 
lattice axes. The larger of the two is called the span- 
ning length of the loop and the smaller one is called its 
roughness. Consequently, a loop is called percolating if 
its spanning length is equal to the system size L. This is a 
binary decision for each disorder instance that is further 
used to obtain the percolation probability Pl(p,Pi/ii) for 
a lattice of a certain size at a given value of the disor- 
der and dilution parameters p and p — Pim, respectively. 
Note that system spanning loops are a direct consequence 
of the periodic boundary conditions we employ. For com- 
pleteness, we note that we used the number 1 to signify 
a sample containing a percolating loop and otherwise. 
In practice, we fix one of the parameters, say p, and tune 
the value of p so as to determine the critical point p c be- 
yond which percolating loops emerge in the limit of large 
system sizes at that particular value of p. For clarity, 
Pl{p,p) at fixed p is simply denoted as Pl{p)- Refer- 
ring to percolation theory j5|, Pl(p) is expected to scale 
as Phip) ~ fi[(p — Pc)L 1 ' 1/ ], where v denotes the criti- 
cal exponent that describes the divergence of the typical 
length scale as the critical point is approached. Since 
we consider string-like objects, it is possible to find more 
than just one percolating loop per disorder instance. In a 
preceding work fl6| we found that the average number of 
percolating loops per sample satisfies the scaling relation 
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FIG. 3: Critical lines in the disorder-dilution planes. Phase 
diagram for type I and type II dilution, (a) and (b) respec- 
tively. (A) signifies a domain in the disorder-dilution plane 
without percolating loops and (B) labels a domain where per- 
colating loops are present. 



(N) ~ } 2 [{p-Pc)L 1/v ) governed by the same values for 
p c and v as the percolation probability. Therein (N) was 
found to exceed a value of 1 above the critical point, sim- 
ilar to what was found earlier for a different loo p p ercola- 
tion model subject to optimization constraints [17] . This 
is in contrast to usual random-bond percolation, where 
the objects of interest are cluster of connected sites on 
the lattice and where there is a unique percolating cluster 
above the percolation threshold. 



A. Type I dilution 

As described above, we allow for a fraction p\ of edges 
that carry zero weight. A loop can include those edges 
without increasing its weight or the configurational en- 
ergy. As a consequence, for values of p\ > percolating 
loops can emerge even if the disorder parameter p has 



TABLE I: Critical points in the disorder-dilution plane for 
type I dilution. From left to right: numerical value of the 
disorder (p) parameter and the dilution (pi) along the criti- 
cal line (the values without errorbars where kept fixed) and 
critical exponent v that describes the divergence of the corre- 
lation length. S measures the quality of the data collapse as 
the mean-square distance of the data sets to the master curve 
in units of the standard error. 
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FIG. 4: Scaling analysis to probe points along the critical 
lines in the disorder-dilution planes, (a) For a disorder pa- 
rameter fixed to p — 0.25, the scaling analysis of the perco- 
lation probability Pl yields a critical value p\ — 0.184(2) for 
the dilution above which percolating loops appear in the sys- 
tem. It further yields a critical exponent v = 1.52(11). (b) 
For a disorder parameter fixed to p = 1.0, the scaling anal- 
ysis of the normalized average box size Al yields a critical 
value Pij = 0.4998(5) for the dilution parameter and a critical 
exponent v = 1.33(5). 



a value smaller than its critical value along the pi = 
axis, i.e. p — 0.34. In particular for p\ > 0.5, there is 
a system spanning cluster of edges with zero weight and 
for L — > oo an arbitrarily small value of p will suffice to 
generate percolating loops on the lattice. 

First, we probed the critical line that separates the 
nonpercolating (small p\ at fixed p) and the percolat- 
ing domains in the disorder-dilution plane. For this 
purpose, we performed simulations for systems of size 
L = 24 . . . 128, the number of samples ranging from 32000 
for L = 24 to 12800 for the largest systems. Within the 
scaling analysis, we fixed the disorder parameter p at val- 
ues listed in table [I] and varied the dilution p\ so as to 
determine the critical point p\ at which percolating loops 
emerge. As an indicator for the transition, we therefore 
study the finite-size-percolation probability P^(pi). The 
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resulting phase diagram is depicted in figure [Ha) an d 
a scaling analysis for one of the points along the criti- 
cal line is shown in figure [4|a). A data collapse involv- 
ing data sets for different system sizes L can be used to 
estimate p\ and v. The quality of the data collapse is 
therein measured by the mean square distance S of the 
data sets to the master curve in units of the standard er- 
ror [l(| , see table |TJ For the example illustrated in figure 
IDJa), where the disorder parameter was fixed to p = 0.25, 
the scaling analysis resulted in a value p c — 0.184(2) for 
the dilution above which percolating loops are present in 
the system. It further yields a value v = 1.52(11) for 
the critical exponent that describes the divergence of the 
correlation length. For consistency we checked some of 
the points by performing simulations along both axis: in 
addition to the data obtained at fixed p we performed 
simulations at fixed p\, probing the critical point along 
the p-axis. For example, probing along the pi-axis at 
p = 0.2 we found p\ = 0.264(2), while probing along 
the p-axis at p\ = 0.264 we obtained p = 0.200(2) in 
complete agreement. Albeit the sample numbers in the 
present study are somewhat smaller than in [l6| that led 
to the estimate of v at p\ — 0, the numerical values we 
obtain for v along the critical line are consistent with the 
value of v = 1.49(7) reported there. This leads us to 
conclude that type I dilution does not change the uni- 
versality class of NWP. Note that it is not possible to 
find negative weighted loops along the p\ = 1 and p = 
axes since at least one edge with a negative edge weight 
is need in order to form a loop. 



B. Type II dilution 

As detailed above, we allow for a fraction pn of absent 
edges on the lattice. In terms of random-bond perco- 
lation, one can expect that the lattice decomposes into 
independent clusters of connected sites for p\\ > 1/2. 
Therefore its very unlikely to find negative weighted loops 



TABLE II: Critical points in the disorder-dilution plane for 
type II dilution. From left to right: numerical value of the 
disorder (p) and dilution (pn) parameter along the critical line 
(the values without errorbars where kept fixed) and critical 
exponent v that describes the divergence of the correlation 
length. 5* measures the quality of the data collapse as the 
mean square distance of the data sets to the master curve in 
units of the standard error. 
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FIG. 5: FSS of the average normalized loop length P£° that 
leads to the estimate of j3 listed in table [TTTl Here, the data 
collapse is best near the critical point while there are correc- 
tions to scaling off criticality. 
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FIG. 6: Analysis of the fractal dimension df, the roughness 
exponent d r and the exponent j/v that relates to the fluctu- 
ations of the loop lengths at p = 1.0 and pn = 0.4998(5). (a) 
Estimation of the above exponents by means of a pure fit to 
the power law data, (b) Estimation of the asymptotic expo- 
nents by extrapolation using the effective (local) exponents of 
the data points shown in (a). 



with a spanning length equal to L for values of the di- 
lution that exceed pu = 1/2. Proceeding as above, we 
fixed the value of the disorder parameter to the values 
listed in table [Til an d tuned the value of the dilution pa- 
rameter so as to probe the corresponding critical points 
Pjj, above which percolating loops cease to exist (large 
pn at fixed p). We performed simulations for systems 
of size L = 32 . . . 128. The number of samples ranged 
from 19200 for L = 32 to 12800 for the largest sys- 
tems. We performed in each case a FSS analysis. Only 
at p = 1.0 we carried out further simulations at L = 192 
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(9600 samples) so as to decrease finite-size effects fur- 
thermore and to establish a more reliable scaling analy- 
sis to settle a value for the critical exponent v. Figure 
H^b) illustrates a scaling analysis for fixed p — 1.0, cor- 
responding to the rightmost data point in figure [3fb). 
Instead of the percolation probability shown for the type 
I dilution in figure dja), we here show the average area 
(Al) = (£ x x £ y )/L 2 of the smallest box that fits the 
largest loop on the lattice normalized by the system size. 
In the same manner as the percolation probability, this 
quantity approaches a step function in the limit of large 
system sizes and its finite size scaling behavior follows 
the analogous form (Al) ~ fsKpu — p^L 1 '"]. For fixed 
p = 1.0 the FSS of the observable Al yields a critical 
value pfi = 0.4998(5) of the dilution parameter and the 
critical exponent v = 1.33(5) with a quality S — 1.25 of 
the data collapse. These findings are in complete agree- 
ment with the corresponding values listed in table HT1 that 
where obtained by means of the percolation probability. 
It is interesting to note that the exponent v steadily de- 
creases from the value v = 1.49(7) found for the undi- 
luted model at p = 0.340(1) to the value v = 1.33(5) 
at p = 1.0 and p u = 0.4998(5), see table HTT1 The lat- 
ter two values are also characteristic for random bond 
percolation in 2d, that is the accompanying process for 
NWP subject to the type of dilution under investiga- 
tion. The change of the exponent v indicates the intro- 
duction of type II dilution changes the universality class 
of NWP. In particular, since the value of v is compati- 
ble with the standard percolation exponent, we ask the 
question whether the NWP is here equivalent to stan- 
dard percolation, although we look for loops of negative 
weight instead for percolating clusters. In the following 
we present a more detailed analysis of the critical point at 
p = 1.0 and pu = 0.4998(5), which shows that indeed the 
NWP behavior at p — 1.0 is not equivalent to standard 
percolation. 



C. Detailed analysis at p = 1.0 and pfj = 0.4998(5) 

Above we found the critical point p n and the exponent 
v from a FSS analysis of the average normalized box size 
Al- Further, the average normalized loop length P£° = 
(£)/L 2 can be used to obtain a further critical exponent, 
namely the percolation strength /?, by use of the scaling 
assumption P£° ~ L~P/ V fi\(pu — pf^L 1 ^] as suggested 
by the corresponding observables known from percolation 
theory. From the FSS of P£° we obtained the critical 
exponent [3 = 0.89(6) with a quality S — 0.97 of the 
scaling assumption, see figure 151 table IIII1 In the formula 
above, the value of I refers to the length of the largest 
loop found for a given realization of the disorder. At 
criticality this reduces to P£° ~ L~P/ U . 

Next, we consider the fractal (scaling) dimension of the 
loops defined via (£) ~ L d f . The value of df is bounded 
by the intrinsic dimension of a line and the topological 
dimension of the environment to 1 < df < 2. From seal- 
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FIG. 7: Detailed analysis at p = 1.0 and pu. Rescaled distri- 
bution Dl(£) of the length I of the largest loop found for each 
realization of the disorder. The inset shows the distribution 
n a of all loop-lengths s, excluding the spanning loops, at the 
critical point. Grey data points were omitted from the fit that 
led to the estimate of r = 2.54(5). 

ing theory, we obtain the relation df = 2 — f3/v. This 
scaling relation may further be used to check for consis- 
tency as we employ two different methods to estimate 
critical exponents, i.e. a data collapse involving data at 
different values of pu and an analysis of the algebraic 
scaling of different observables right at the critical point. 
A simple fit to the power law data, shown in figure Oa), 
yields the exponent df = 1.333(2) with quality Q — 0.51 
8]. This is compatible with the scaling relation within 
error bars, given the values deter minded above for (3 
and v. Interestingly, the value of df found here coincides 
with the fractal dimension of self avoiding walks on a 
regular lattice (d^ AW = 4/3) and of exterior perimeters 
measured for 2d percolation clusters at the percolation 
threshold Moreover, it compares well with the value 
dJ B = 1.34(2) that describes the scaling of the average 
backbone perimeter with system size, measured for ran- 
dom percolation at the percolation threshold 18]. This 
suggests, in analogy with the conclusions of [l8|, that 
negative weighted loops at (p,pu) = (1.0,0.4998) belong 
to the same universality class as self-avoiding walks. The 
roughness of the loops is expected to display the critical 
scaling (r) ~ L dr with a roughness exponent d r = 1 . The 
fit to the power law data, also shown in figure^a), yields 
indeed the exponent d r = 0.998(2) with quality Q = 0.33, 
highlighting the self-similar nature of the loops. 

The susceptibility like quantity xl = ((£ 2 ) - (£) 2 )/L 2 , 
that measures the mean square fluctuation of the loop 
length, exhibits the critical scaling \l ~ L 1 ^ that allows 
to determine the critical exponent 7 = 0.89(4) with a 
quality Q — 0.56 of the fit, see again figure [U(a). It is 
connected to the other exponents via the scaling relation 
7 + 2/3 — 2v, which is fulfilled within error bars. 

Note that an analysis of the effective (local) exponents 
of the data belonging to (£) , (r) and \l yields the asymp- 
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TABLE III: Critical points and exponents. From left to right: Type (RBP: random bond percolation, NWP: negative weight 
percolation ), dilution (pn) and disorder (p) parameters (note that there is no entry for RBP since the disorder parameter p 
introduced here has no analog in usual random percolation), critical exponent of the correlation length v, percolation strength 
/3, exponent 7, Fisher exponent r and fractal dimension df at criticality. 



Type 


pn 


Pc 


V 


P 


7 


r 


df 


RBP 


0.5 




4/3 


5/36 


43/18 


187/91 


91/48 


NWP 


0.0 


0.340(1) 


1.49(7) 


1.07(6) 


0.77(7) 


2.59(3) 


1.266(2) 


NWP 


0.4998(5) 


1.0 


1.33(5) 


0.89(4) 


0.89(6) 


2.54(5) 


1.333(2) 



totic scaling exponents df = 1.33(1), d r = 0.99(1) and 
7 = 0.89(6) compatible with the values found above but 
with somewhat larger errorbars, see figure [H^b). 

As evident from figure [71 the distributions Dl(£) of the 
largest loops £ found for each realization of the disorder 
for different system sizes L show a nice data collapse 
after a rescaling to the form D L {1) = L- d f f 5 [£/L d f}. 
This means that not only the average {£) but the full 
distribution scales with df, again highlighting the fractal 
nature of the loops. Therein, the peak close to £/L d f = 
0.25 stems from the non-spanning loops and the peak 
close to £/L df =0.75 stems from those that truly span 
the lattice along at least one direction. At the critical 
point, we further expect the distribution n s of all loop 
lengths s, excluding loops with a spanning length L, to 
exhibit an algebraic scaling n s ~ s~ T , where r signifies 
the Fisher exponent. The exponent r characterizes the 
complete ensemble of loops and it is related to the fractal 
exponent of the loops via r — 1 = d/df. Here, we found 
a value of r = 2.54(5), see inset of figure [Jj in agreement 
with the latter scaling relation. 

IV. CONCLUSIONS 

In summary, we performed numerical simulations in 
order to probe the effect of dilution on the critical prop- 
erties of NWP on 2d lattice graphs. In this context and 
by means of exact combinatorial optimization algorithms 
we have investigated two different types of dilution. One 
type affecting the distribution of weights that are as- 
sociated with the edges, the other type having an im- 
pact on the topology of the underlying lattice. We used 
observables from percolation theory and a FSS analysis 
so as to probe critical points along the critical lines in 



the disorder-dilution plane that separate domains that 
allow/disallow system spanning, i.e. percolating, loops. 
We found that the first type of dilution does not change 
the universality class of NWP, i.e. the critical exponents 
are the same as for the undiluted model. The second 
type of dilution very well leads to a change of the uni- 
versality class: As the disorder on the lattice increases, 
the critical exponent v as well as the critical point of 
dilution approach the corresponding values known from 
usual random bond percolation, although the process, 
also due to its string-like nature, is clearly different from 
standard percolation, also visible through the values of 
the other critical exponents. A more detailed analysis 
of the critical exponents for one particular critical point 
in the disorder-dilution plane, i.e. (p,pn) = (1.0, 0.4998), 
verified that the exponents are connected by the usual 
scaling relations. Moreover, the results for the scaling 
behavior of the percolating loops suggest that negative 
weighted loops at that particular critical point belong to 
the same universality class as self-avoiding walks. 
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